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Abstract 

A search is described for the neutral Higgs bosons h° and A*^ predicted by models with two 
scalar field doublets and, in particular, the Minimal Supersymmetric Standard Model (MSSM). 
The search in the Z°h° and h^A" production channels is based on data corresponding to an 
integrated luminosity of 25 pb~^ from e"'"e~ collisions at centre-of-mass energies between 130 and 
172 GeV collected with the OPAL detector at LEP. The observation of a number of candidates 
consistent with Standard Model background expectations is used in combination with earlier 
results from data collected at the Z*^ resonance to set limits on mh and in general models 
5^ , with two scalar field doublets and in the MSSM. For example, in the MSSM, for tan/5 > 1, 
minimal and maximal scalar top quark mixing and soft SUSY-breaking masses of 1 TeV, the 
95% confidence level limits mh > 59.0 GeV and ttia > 59.5 GeV are obtained. For the first 
time, the MSSM parameter space is explored in a detailed scan. 
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1 Introduction 



The data accumulated with the OPAL detector at centre-of-mass energies of 130 — 172 GeV, 
corresponding to approximately 25 pb~^ of integrated luminosity, have opened up new kinematic 
domains for particle searches. We describe searches performed in these high energy data for 
neutral Higgs bosons which decay to hadrons or tau-leptons. 

In the Standard Model (SM) |l[], spontaneous symmetry-breaking is effected by the self- 
interaction of one scalar (Higgs) field doublet [Q. The model predicts one Higgs boson, Hgj^, 
the mass of which is not specified. The OPAL search for B.^-^ at ^/s = 130 — 172 GeV has 
resulted in a lower bound on its mass of ^h°^ > 69.4 GeV at the 95% confidence level (CL) 0. 

Supersymmetry (SUSY) |Q is considered an attractive possible extension of the SM, since 
it provides a solution to one of the outstanding problems of the SM, that of the "hierarchy" of 
energy scales ^ . The implementation of SUSY requires at least two Higgs field doublets. There 
are exactly two in the Minimal Supersymmetric extension of the Standard Model (MSSM) 0. 
These fields couple separately to up- type quarks for the first doublet, and to down-type quarks 
and charged leptons for the second doublet, and have vacuum expectation values vi and V2, 
respectively. Scalar field doublets that couple in this manner may exist more generally and the 
class of such models is known as Type II Two Higgs Doublet Model (2HDM) 0. Here, 2HDM 
models will be understood to have no extra particles besides those of the SM and the two 
scalar doublets, and there is no mass relation between the different neutral and charged scalar 
field particles. In any two Higgs field doublet model, the Higgs sector comprises five physical 
Higgs bosons: two neutral CP-even scalars h'' and H° (with masses satisfying < m-n by 
definition), one CP-odd scalar A*^ and two charged scalars H^. In this paper the searches are 
restricted to the neutral Higgs bosons h*^ and A°. The heaviest neutral Higgs boson H° is likely 
to have a mass beyond the reach of LEP, and OPAL searches for H^ bosons have been published 
separately [§]. 

At the current e"'"e~ centre-of-mass energies (a/s) accessed by LEP, the h° and A'' bosons 
are expected to be produced predominantly via two processes: the "Higgs-strahlung" process 
e"'"e~— >-h°Z° (where the Z° boson is on-shell) and the "pair production" process e"'"e~— > 
Contributions from the W"'"W~ and Z°Z° fusion processes account for a small part of the total 
production, except close to the kinematic limit of the e^e~— ^h'^Z" process 0. For these two 
principal processes, the cross-sections a^z and cThA are related at tree-level to the SM cross- 
sections [ pil[| : 

e+e--.hOzo : atz = sm\p - a) (1) 
e+e--^hOAO : = cos\(3 - a) X a^^, (2) 

where and are the cross-sections for the SM processes e"'"e^^Hgj^Z° and e"'"e~— >-i/z/, 
and A is a kinematic factor, depending on mh, m\ and ^/s, typically having values between 0.5 
and 0.7 for the centre-of-mass energies under consideration. The angle jS is defined in terms 
of the vacuum expectation values of the two scalar fields, tan/? = V2/V1, and a is the mixing 
angle of the CP-even (h°,H'') fields. The coefficients sin^(/3 — a) and cos^(/9 — a) indicate 
complementarity in the cross- sect ions of the two processes, a feature which is exploited in 
deriving bounds for Higgs boson masses and other model parameters. 

In the MSSM at tree level the following mass relations are predicted: mh < mz, < ^h, 
< '^H and mH± > mw± 0- Loop corrections, dominantly from the top and scalar top 
quarks (t), strongly modify these mass relations and also have some moderate impact on the 
Higgs boson couplings |jll|. The shift in mh, approximately proportional to m^ x log(m5:^/m^). 
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can be several tens of GeV. While the top quark mass has been measured to be mt 



(175±5) GeV |T^, the mass of the scalar top quark mi depends on the mixing in the t sec- 
tor which, in turn, depends on several other parameters of the MSSM. However, even for a 
choice of SUSY parameters which maximises the mass shift, mh should not exceed approxi- 
mately 135 GeV |Ty]. Although this bound is beyond the ultimate reach of the LEP collider, a 
substantial fraction of this mass range is accessible. 

In this work we undertake a more detailed examination of the MSSM parameter space 
than has been done in the past. The coefficients sin^(/? — a) and cos^(/5 — a) which appear in 
Eqs. (|ip and (|^) depend on a number of MSSM parameters which enter via mixing in the t 
sector and which will be summarised in a later section. We perform detailed scans over broad 
ranges of these parameters. Each of these scans is considered as an independent "model" within 
the MSSM, and results are provided for each. In increasing order of generality, they include: 
(A) Particular choices of parameters for "minimal" and "maximal" t-mixing and soft SUSY- 
breaking masses fixed to be large as defined in |T^ which are considered as the "benchmark 
case" and adopted by most search groups, (B) varying the parameters over a wide range, but 
keeping relations between some of them corresponding to minimal and maximal mixing, and 
(C) a "general" scan where the MSSM parameters are allowed to vary independently within 
wide, but reasonable, ranges. 

The final-state topologies of the processes (0) and (|^) are determined by the decays of the 
Z°, h° and A° bosons. Since the Higgs bosons couple to fermions with a strength proportional 
to the fermion mass, the Higgs dominantly decays into pairs of the most massive fermions 
which are allowed by the kinematics, most notably b quarks and tau-leptons for the LEP mass 
range. For particular choices of the model parameters (e.g. for tan (3 < 1) decays into cc 
may be enhanced. For 2mA < '^h, the process h°— i>A°A'' is also allowed and may even be 
the dominant decay, leading to more complex final states than those from direct decays into 
fermion pairs. In the MSSM, Higgs bosons may also decay into SUSY particles if allowed by 
kinematics. In particular, the decay into pairs of neutralinos may lead to "invisible" Higgs 
decay channel^ which must be considered in a full treatment of the model. 

In searching for the process e"'"e~— >-h°Z°, the fact is exploited that in most of the MSSM 
parameter space with tan f3 > 0.7 the decay properties of the h*^ boson are essentially those of the 
SM Higgs boson. Thus, earlier OPAL searches for the Hg^^ boson 0, including those performed 
at energies above the Z° mass [|13|, [1^, are interpreted here as searches for e"'"e~— >h°Z°. The 
reduction of the search sensitivity due to sin^(/9 — a) in Eq. (1) is taken into account. Dedicated 
searches for "invisible" final states at a/s ~ rnz [0 are also included. Since these searches are 
published, we only summarise the results which are relevant for the present purpose. These 
searches for Hg^^ are also efficient for e+e~-^h°Z°^A°A°Z°, sometimes after small modifications 
to the search. 

In searching for the process e'^e~^h.^A^, the following final states are most important: 
(h*^— i>bb)(A°— >-bb), (h°— >r"'"r~)(A''— i>qq) and (h"— >-qq)(A°— i>r"'"r~). The searches in these chan- 
nels using the data accumulated by OPAL between ^/s = 130 and 172 GeV have not been pub- 
lished. Therefore they are described in greater detail. The search for h°A'^^A'^A°A°— s>bbbbbb 
is important when the decay h^^A^A^ is kinematically allowed and is also presented here for 
the first time for data taken above the Z" energy. 

This paper starts in Section ^ with a short description of the OPAL detector, the data 
samples used and the various Monte Carlo simulations used to obtain the detection efficiencies 

^In the MSSM R-parity is conserved and throughout this paper we assume that the hghtest supersymmetric 
particle is the Hghtest ncutrahno, x^. 
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and to estimate the backgrounds from SM processes. This is followed by a description of the 
event selections for the various h°Z° and h°A° channels in Sections ^ and 

A new statistical method |T^, summarised in Section ^ has been used to combine the results 
from different search channels and data sets taken at different centre-of-mass energies. 

The model-independent and 2HDM results are summarised in Section ^, followed by those 
interpreted within the MSSM in Section |^. 

Previous OPAL searches for the h° and A° bosons, based on data collected at y/s^mz are 
described in [pTSI , [l^. The relevant publications from the other LEP collaborations describing 
neutral SUSY Higgs boson searches are listed in ]T7| . 



2 Experimental considerations 

The present search includes data collected with the OPAL detector in 1995 at y/s = 
130 - 136 GeV (5.2 pb~^), in 1996 at = 161 GeV (10.0 pb"^), and at 170 - 172 GeV 
(10.4 pb"^). The resuhs are combined with those from earlier searches [1111 which are based on 
the analysis of up to 4.5 million hadronic decays. 

The OPAL detector is an apparatus with nearly complete solid angle coverage and excellent 
hermeticity. The central tracking detector consists of a high-resolution silicon microstrip vertex 
detector (/iVTX) [l^ with polar angle| coverage | cos6'| < 0.9, which immediately surrounds the 
beam-pipe, followed by a high-precision vertex drift chamber, a large- volume jet chamber, and z- 
chambers, all in a uniform 0.435 T axial magnetic field. A lead-glass electromagnetic calorimeter 
is located outside the magnet coil, which, in combination with the forward calorimeter, gamma 
catcher, and silicon-tungsten luminometer pOf, complete the geometrical acceptance down to 



33 mrad from the beam direction. The silicon-tungsten luminometer serves to measure the 



integrated luminosity using small-angle Bhabha scattering events |^ . The magnet return yoke 
is instrumented with streamer tubes for hadron calorimetry and is surrounded by several layers 
of muon chambers. 

Events are reconstructed from charged-particle tracks and energy deposits ("clusters") in 
the electromagnetic and hadronic calorimeters. The tracks and clusters must pass a set of 



quality requirements similar to those used in previous OPAL Higgs boson searches p2|. In 
calculating the total visible energies and momenta, i^vis and Pvis? of events and of individual jets, 
corrections are applied against double- counting of energy in the case of tracks and associated 
clusters ||2^, |2^. For the analysis presented here, charged particles and neutral clusters are 
grouped into jets using the Durham algorithm ^5 . 



The tagging of jets originating from b quarks plays an important role in Higgs boson 
searches, since both h'' and A'' decay preferentially to bb over large domains of the two-field- 
doublet and MSSM parameter spaces. Primary and secondary vertices are reconstructed in 
three dimensions following two algorithms, using only tracks which pass an additional set of 
quality requirements. The first method (BTAGl) considers all such tracks in a jet and 
attempts to fit them to a common vertex. Tracks are discarded from the vertex by an iterative 
procedure which drops the track with the largest contribution to the vertex fit, until the 
largest contribution is less than 4, with at least two tracks remaining. In the second method 
(BTAG2) [^, the intersection of all pairs of such tracks in a jet having impact parameter 

^ OPAL uses a right-handed coordinate system where the +z direction is along the electron beam and where 
+x points to the centre of the LEP ring. The polar angle, 6, is defined with respect to the +z direction and 
the azimuthal angle, 0, with respect to the horizontal, +x direction. 
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significance h/cjh (where h is the impact parameter^ and its error) greater than 2.5 are con- 
sidered as seed vertices. The other tracks in the jet are added to the seed vertex in the order 
which results in the greatest vertex probabihty after each addition. The process continues until 
either all tracks in the jet are added or the resulting vertex probability falls below 1%, in which 
case the last track to have been added is dropped. If more than one acceptable vertex per jet 
is found via this algorithm, the best is chosen according to a set of criteria |26| involving track 
multiplicity and the vertex decay length significance S = L/a^, where L is the vertex decay 
lengthQ and ai its error. In both algorithms, S is then used to distinguish between b-flavoured 
hadron decays and background. The methods are found to not be fully correlated, and using 
them both adds discriminating power. Charged track impact parameters are also used to com- 
plement the secondary vertex algorithms via the forward multiplicity, defined as the number 
of tracks in a jet with b/ab > 2.5. Finally, semi-leptonic b-hadron decays are exploited by first 
identifying electrons via the method described in ||28| and muons with the algorithm described 
in and then considering the transverse momentum with respect to the corresponding jet 
axis. 

The signal detection efficiencies and accepted background cross-sections are estimated using 



a variety of Monte Carlo samples. The HZHA generator |^ is used to simulate Higgs boson 
production processes. The detection efficiencies are determined at fixed values of the Higgs 
boson masses using sample sizes varying between 500 and 10,000 events. Efficiencies at arbi- 
trary masses are evaluated using spline fits in the (mh,mA) plane between these points. The 



background processes are simulated primarily by the following event generators: PYTHIA ^1 
((Z77)*^qq(7)), EXCALIBUR H and grc4f [||] (four-fermion processes (4f)), BHWIDE ^ 



(e+e-(7)), KORALZ |35| (/i+/i"(7) and r+r-(7)), PHOJET |3| and Vermaseren (hadronic 
and leptonic two-photon processes (77)). The generated partons are hadronised using JET- 
SET 1^ and the resulting particles are processed through a full simulation |^ of the OPAL 
detector. 



3 Searches for the process e+e — ^h Z 

The OPAL searches for the SM process e"'"e~— >IIgi^Z*^ are interpreted as searches for the process 
e'^e~—>-h^Z^ via Eq. (1). The searches for Hgj^ using all data recorded at center-of-mass energies 
up to 172 GeV are described in [Q. They make use of the following final states: 

. "Four-jet": (HOM-^bb)(ZO^qq) (q=u,d,s,c,b), 

• "Missing energy" : (IIgi^-^qq)(Z°— >z/z/) (q includes quarks and gluons), 

• "Charged lepton": (Hgj^— i>qq)(Z°^e+e^, /i+yU^), and 

• "Tau-lepton": (HgM^qq)(Z°^r+r-) and {Rl^^T+T-){Z^^qq). 

The search in the missing energy channel is also sensitive to small contributions to H^m 
production coming from the W"'"W~ fusion process e"'"e~^z/z/Hgf^ while the search in the 

The impact parameter is taken to be positive if in the two-dimensional projection the track path crosses 
the jet axis in the direction of the flight direction; otherwise it is negative. 

^ The vertex decay length is the projection onto the jet direction of the distance between the primary vertex, 
as reconstructed for the event (see for the algorithm), and the secondary vertex, as reconstructed for the 
jet. The decay length is taken to be positive if the vector that connects the primary to the secondary vertex is 
at an angle of less than 90° from the direction of the associated jet, and negative otherwise. 
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15 


19 
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Table 1: Summary of the searches for the SM Higgs boson at centre-of-mass energies of 161 and 
170 — 172 GeV. For each channel the signal detection efficiencies for ^ng^ = 60 and 70 GeV, 
the number of expected background events, and the number of events selected are given. The 
statistical error on the efficiencies is typically 1-4%. 



charged lepton channel is sensitive to those from the Z^Z" fusion process e"'"e~— *-e"'"e~Hg]yj. 
These contributions are taken into account in the corresponding channels. 

The results from these published searches at = 161 — 172 GeV are summarised in 
Table |I], which lists the signal detection efficiencies for two Higgs boson masses in the range 
of interest, the residual expected number of background events and the number of selected 
candidate events in each search channel. As can be seen from the table, the selection criteria 
applied at ^/s = 161 GeV select one candidate in the missing energy channel, with a mass 
of (39.3±4.9) GeV, while those applied at 170 — 172 GeV select one candidate in the four- 
jet channel with a mass of (75.6±3.0) GeV. The earlier OPAL searches applied to Z^ boson 
decays [0 selected one candidate in the charged lepton channel /i+/i~Hgp^ with a mass of 
(61.2±1.0) GeV (with 0.38±0.04 events expected from background) and two candidates in the 
missing energy channel with masses of (6.3±0.8) GeV and (24.8±3.0) GeV (with 2.3±0.4 events 
expected from background) . All these candidates are considered as possible Higgs boson events 
when limits in the MSSM parameter space are computed. 

The above searches are also sensitive to the process e'^e~^h^Z^ followed by h°— i>A°A°. The 
selection is slightly changed with respect to that described in in the case of the four-jet 
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Table 2: Signal detection efEciencies for the searches for the SM Higgs boson, at centre-of-mass 
energies of 161 and 170 — 172 GeV, apphed to the processes with h*^— >A°A° followed by A°^bb. 
The efficiencies are quoted for mh = 60 GeV and rriA = 30 GeV, with typical statistical errors 
of 1-4%. 



and missing energy channels. In the four-jet channel the likelihood discriminant, as described 
later in Section ^TT], is reoptimised for the Z^A^A*^ final state. This increases the efficiency from 



23% to 29%, while the background expectation remains approximately equal. In the case of 
the missing energy channel at ^/s = 161 GeV a requirement on the jet resolution parameter, 
?/23 < 0.05, limits the acceptance to two-jet events only. Removing this requirement, the 
detection efficiency for h°— i>A°A° (mh = 60 GeV, ttia = 30 GeV) increases from 38% to 68% 
while the background increases from 0.9 to 1.1 events. For the decay of the A° boson, only 
the predominant bb final state is considered. For the charged lepton and tau-lepton channels, 
Monte Carlo simulations have demonstrated that the detection efficiencies for this two-stage 
process are close to those of the h° decay to fermion pairs. The detection efficiencies for the 
particular case of mh = 60 GeV and = 30 GeV, a point close to the boundary of the 
kinematically- allowed region for h°— >A°A°, are shown in Table |[ 

4 Searches for the process e+e~^h^A^ 

In this section the searches for the MSSM process e+e~^h'^A'^ for final states (h°-^bb)(A'^— >-bb), 
(h°^r+r")(A°^qq), (h°^qq)(A°^r+r") and (h°^A°A°^bbbb)(A°^bb) are summarised. 
If the same branching ratios to bb and r+r^ are assumed for h° and A° as for Hgf^, then the first 
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three of the above final states account for approximately 90% of all h°A° decays, in the portion 
of the {m\,mii) phase space where > frii^. For points in this phase space where mh > 
and the decay h°— i>A°A° is kinematically allowed, the last final state represents approximately 
66% of decays of this type. 



4.1 The channel h^A^^bbbb 

The signature for events from the process h°A°— i>bbbb is four energetic jets containing b- 
hadrons and a visible energy close to the centre-of-mass energy. As ^/s changes from 130 GeV 
to 172 GeV, the background changes considerably in size and composition. At 130 GeV and 
161 GeV the main background comes from (Z°/7)*— >qq with or without initial state radiation 
accompanied by hard gluon emission. Four-fermion processes, in particular e"'"e~^W"'"W~, 
play a minor role since the threshold for these processes is at most only marginally crossed. 
At 170 — 172 GeV this background becomes more important, while the (Z°/7)* background is 
reduced. Two-photon processes have a large cross-section at all energies; however, the selection 
requiring multihadronic final states with a visible energy close to ^/s reduces them to a negligible 
level. 

The selection proceeds in two phases. First, a preselection is applied to retain only those 
events which have some similarity to the signal. The events remaining after preselection are 
then analysed using a likelihood technique. 

The preselection consists of the following requirements: 

(1) The events must qualify as being hadronic final states as described in [ ^ . 

(2) The radiative process e"^e~— >(Z°/7)*^qq7 is largely eliminated by requiring that the 
effective centre-of-mass energy, y/s', obtained by discarding the radiative photon from 
the event following |0|, is greater than 110, 140 and 150 GeV for ^ = 130 - 136, 161, 
and 170 — 172 GeV, respectively. 



(3) The events are reconstructed into four jets using the Durham algorithm ||2^. The jet 
resolution parameter 1/34, at which the number of jets changes from 3 to 4, is required to 
be larger than 0.005. The distribution of log 7/34 is shown in Fig. |l|(a). 

(4) The value of the event shape C-parameter must be greater than 0.45. The C-parameter 
is defined as C = 3(eie2 + 6263 + 6361), with ei, 62 and 63, being the eigenvalues of 



the normalised momentum tensor of the event |41|. It ranges from C = for perfectly 



back-to-back two-jet events to C = 1 for perfectly spherical events. 

(5) Each of the four reconstructed jets must contain at least two charged tracks and at least 
two electromagnetic calorimeter clusters. 

(6) A fit of the jet four-momenta, in which the energy and momentum of the final state is 
constrained to that of the initial e"^e~ state, must yield a ^^-probability larger than 0.01. 

The results of the preselection are listed in Table |^. Except for the number of events retained 
after (1), the agreement between observed events and expected background predicted by the 
Monte Carlo simulation is good. The discrepancy after the first requirement can be explained 
by inaccurate modelling of events that radiatively return to the Z° and by not including the 
two-photon events in the Monte Carlo prediction. After a cut on \/s', these events are rejected 
and the background prediction from Monte Carlo describes the data well. 
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Next, a likelihood technique is used to classify the remaining events as either (Z°/7)*— i>qq(7) 
(1), a four-fermion process (2), or A^h"— »bbbb (3). Three kinematic and six b-tagging variables 
are input to the likelihood. The kinematic quantities are: the smallest angle between any jet 
pair, the logarithm of the probability for a fit of the jet momenta with energy, momentum and 
equal di-jet mass constraint, and the smallest di-jet mass difference after the energy momentum 
conserving fit. For the b-tagging quantities used here, all vertices on which their calculation is 
based must contain at least two tracks, each having two hits in the r — and r — z ladders of the 
/xVTX detector. The calculation of forward multiplicities takes into account only tracks which 
have two hits in the r — (p and r — z ladders of the /iVTX detector. These quantities are (see 
Section ^ for a description) : the sum of the two largest and the sum of the two smallest decay 
length significances in the four jets for vertices reconstructed with algorithm BTAGl, the sum 
of all decay length significances for vertices reconstructed with algorithm BTAG2, the sum of 
the two largest and the sum of the two smallest forward multiplicities, and the sum of the two 
largest transverse momenta with respect to the corresponding jet axis for identified leptons. For 
all of these quantities there is good agreement in their distributions between data and Monte 
Carlo background predictions, and as examples the distributions for the sum of the two largest 
lepton transverse momenta and the sum of the two largest decay length significances for the 
vertex algorithm BTAGl are shown in Fig. |T](b) and (c), respectively. 

For each of these input variables (labelled by i) a normalised histogram fi{xi) is constructed 
from the Monte Carlo for each of the two background classes and the signal (labelled by 
j = 1, 2, 3). For a single variable, the probability for an event to belong to class j is: 

Piixi) = ^ 



and the joint discriminating variable for class j is defined as: 

J2k=i Y\.i=i Pi i^i) 

where the product runs over the nine input variables. The signal likelihood is defined as: 



and is required to be greater than 0.8 for the final selection. The likelihood distribution is 
shown in Fig. [I|(d) and over the entire range good agreement between the data and Monte 
Carlo background prediction can be observed. Some irreducible background from four-fermion 
processes also shows up as a small peak near a likelihood of 1. In total one event is selected, 
with a centre-of-mass energy of 172 GeV, while 1.65, 1.61 and 1.95 events are expected from 
the background simulations at ^/s = 130 — 136, 161, and 170 — 172 GeV, respectively. The 
Poisson probability to select one or fewer events when 5.2 are expected is 2.3%. 

The signal detection efficiencies are affected by the following main uncertainties, all ex- 
pressed as relative percentages. The error from Monte Carlo statistics is typically 4 — 10%. 
However, near the region where the limit is set larger Monte Carlo samples were generated, 
and a fit is made through the grid of efficiencies in the (mh, plane, so the efficiency at a 
certain mass point is effectively based on higher statistics. This results in a statistical error 
of approximately 1% in the region near the limit. The preselection requirements on y/s', y^^, 
and C were varied by amounts equivalent to the difference between the mean value of data 
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Figure 1: Selection variables relevant for the h°A°— i>bbbb analysis, (a) The logarithm of the 
jet resolution parameter 7/34 after preselection requirements (1) and (2). (b) The sum of the 
two largest lepton transverse momenta after the preselection, (c) The sum of the two largest 
secondary vertex decay length significances of the BTAGl algorithm (defined in Section ^ 
after the preselection, (d) The A°h° likelihood distribution after the preselection. The three 
centre-of-mass energies are added for all histograms. The points represent the data. The 
shaded histograms show the qq background and the open histograms show the four-fermion 
background, normalised to the integrated luminosity of the data. Two-photon processes are 
not included. The dashed line histogram shows the expectation for A^h*^ signal events with 
rriA = mh = 55 GeV, where the displayed production cross-sections have been chosen for 
visibility. 
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(c) 170-172 GeV 
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(b) 161 GeV 
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Figure 2: (a)-(c) The efEciencies (in %) for the detection of h°A°—i>bbbb at the three centre- 
of-mass energies in the (m\^, mpj plane. The numbers are shown only for mp^ > mh due to 
the symmetry in mh and mp^. (d) The position in the (m^, m,pj plane of the six possible mass 
combinations of the candidate event superimposed on the expected Standard Model background 
(scaled up by a factor of 100) for all center-of-mass energies. 
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Mass resolution curves for h^A^— >bbbb. (a) The mass sum of mh and for 
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Figure 3: 
mh + m,A 

of mh and mA for mh — mA = 0. (d) The mass difference for |mh — mA| = 50 GeV. Only the 
combination with the smallest difference between measured and true mass sum or difference is 
plotted. The histograms show the simulated distributions and the solid lines represent smooth 
fitted functions. 
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Table 3: Effect of the selection criteria on data, background (normalised to the integrated 
luminosity of the data) and signal simulation (m^ = niA = 55 GeV) at the three centre-of- 
mass energies for the signal channel h°A°— >bbbb. The quoted errors on the background are 
statistical. 
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and Monte Carlo. This results in an uncertainty estimate for the modelling of the preselection 
variables ranging from 3.5% to 7.6%, depending on the centre-of-mass energy. The uncertainty 
associated with the requirement on ^^-probability was determined to be 2%. To assess the sys- 
tematic uncertainty from the b-tagging, the impact parameter resolution was varied by 10%, 
and the b-hadron charged decay multiplicity was varied by 0.35 units ||^. The resulting errors 
range from 3.6% to 4.5% depending on the detector configuration (data taken in 1995 versus 
data taken in 1996.) The effect of binning in the likelihood was investigated by using a linear 
bin-to-bin interpolation. This gives rise to an uncertainty ranging from 0.0% to 0.4%. The 
theoretical uncertainty on the cross-section is estimated to be 1%. The systematic error on the 
integrated luminosity is 0.6% [^. The total systematic uncertainty is calculated by adding the 



above uncertainties in quadrature. This gives uncertainties ranging from 5% to 9% depending 
on the centre-of-mass energy. 

The residual background estimate has a statistical error of 20%, 7%, 5% at the 130 — 136, 
161, 170 — 172 GeV centre-of-mass energies, respectively. The uncertainty from modelling the 
preselection was estimated in the same manner as for the efficiency, giving uncertainties ranging 
from 3.3% to 16%. The effects of impact parameter resolution and binning in the likelihood 
were also investigated in the same manner as for the efficiency. Uncertainties in the impact 
parameter resolution give errors between 15% and 23%. The effect of modelling the variables 
used in the likelihood selection is cross-checked by reweighting the Monte Carlo events in such 
a manner as to better approximate the data distributions within the ability of the Monte Carlo 
samples to represent the data. It should be noted that the agreement in these distributions 
before the reweighting procedure already gives a good probability. This procedure yields 
errors ranging from 0.4% to 22%, and the corresponding error estimates on the signal detection 
efficiencies are negligible. In all cases errors due to the tracking resolution uncertainty were 
larger than those due to the reweighting procedure, and therefore the former are taken as 
the systematic errors for these two checks. The effect of binning on the likelihood for the 
background was estimated to range from 0.1% to 4.0%. The effect of modelling the Standard 
Model physics was investigated by comparing different Monte Carlo event generators. For the 
Z*^/7*— >qq background the results using the PYTHIA generator were compared to those using 
the HERWIG generator |^3|]. For the four-fermion background the EXCALIBUR and grc4f 



generators were compared. The uncertainties for the physics modelling range from 3.6% to 
14%. The Monte Carlo generators have an uncertainty on the calculated cross-sections of 0.5%. 



Including the error on the integrated luminosity of 0.6% |£T], total relative uncertainties of 
27%, 28%, 24% are assigned to the residual background estimates at center-of-mass energies of 
130 - 136, 161, 170 - 172 GeV, respectively. 

To make use of the mass information in the calculation of exclusion limits, the hypothetical 
Higgs masses have to be determined. The four jets can be combined in three ways into two 
jet-pairs. The invariant masses of all jet-pair combinations of an event are considered. They 
are calculated using a kinematic fit assuming energy and momentum conservation (4C-fit). The 
mass distributions have a non-Gaussian shape. Since h" and A" cannot be distinguished, the 
mass difference AM = |mA — m^l and the sum M = ttia + are considered instead of the 
masses themselves. Figure ^ shows examples of these mass distributions at y/s = 172 GeV. The 
resolution in M is roughly 3 GeV and does not vary with mass. The tails of the distribution of 
M, caused by the reconstructed jets not corresponding to the parton final states, exist largely 
above the mass peak for M = 80 GeV, whereas for M = 130 GeV they are present both above 
and below the nominal mass value. The shape of AM is independent of the mass sum. Its 
resolution is also approximately 3 GeV and the distribution shows substantial tails. Smooth 
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functions fitted to these mass distribution histograms are used for exclusion hmit calculations 
instead of the histograms themselves. 

4.2 The channel h°A°^T+T~qq 

The r^r^ qq final state can be produced via the processes e"'"e~— > -^T^T qq and qqr"''r . 
These processes are characterised by a pair of tau-leptons and a pair of energetic hadronic 
jets. The backgrounds are predominantly from (2^/7)*— >-qq(7) and four-fermion processes. 
The search in these channels was restricted to the data recorded at the centre-of-mass energies 
= 161 and 170 - 172 GeV. 

The selection begins with the identification of tau-leptons, identical to that in [H, using three 
algorithms which classify each tau-lepton candidates as decaying into an electron, a muon, or 
hadrons. 

In the selection that follows, the tau-lepton direction is approximated by that of the visible 
decay products. If two tau-lepton candidates have momentum vectors separated by less than 
23°, one being identified as a leptonic (electron or muon) decay and one as hadronic, the leptonic 
decay is chosen. The following selection, which is identical to that in [Q up to and including 
(4), was made: 

(1) Events are required to have at least two tau-lepton candidates, each with electric charge 
|g| = 1, and at least nine charged tracks. 

(2) Most of the two-photon and ^{7P /'-))* background events are eliminated by requiring 
that the energy in the forward detector, gamma catcher, and silicon-tungsten luminometer 
be less than 4, 10, and 10 GeV, respectively, that | cos6'jniss| < 0.97 and that > 3 GeV, 
where 6'miss is the polar angle of the missing momentum vector and is the total 
transverse momentum of the event. In addition, the scalar sum of all track and cluster 
transverse momenta is required to be larger than 40 GeV. Accelerator-related backgrounds 
in the forward detectors which have not been fully simulated are taken into account via 
small corrections to the signal detection efficiencies. 

(3) The remaining (Z°/7)*— i>qq(7) background is partially suppressed by requiring that events 
contain at least four jets, reconstructed using the cone algorithm as in [||, where sin- 
gle electrons and muons from tau-lepton decays are allowed to be recognised as low- 
multiplicity "jets". Events with an energetic isolated photon are removed, where an 
energetic isolated photon is defined as an electromagnetic cluster with energy larger than 
15 GeV and no track within a cone of 30° half-angle. 

(4) To suppress the process W"'"W~— ^£z/qq', events are rejected if they contain any track or 
cluster with energy exceeding 0.3\/s. Fig. §(a) shows the distribution of the energy of 
the most energetic electromagnetic cluster scaled by ^i, prior to this cut, for the data 
recorded at ^/s = 161 — 172 GeV, the expected backgrounds, and a simulated Higgs boson 
signal with mh = rriA = 55 GeV. 

(5) The three tau-lepton identification algorithms identify 2.3 r candidates per signal event on 
average. Fake candidate pairs are removed by requiring that the sum of the track charges 
be zero and that the candidates satisfy a pairwise isolation requirement, | cosoj-cosajl < 
0.5, where is the angle between the direction of the i-th t candidate and that of the 
nearest track not associated with it. The indices i,j run over all r candidates with i 7^ j. 
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Fig. |^(b) shows the distribution of this variable for pairs of tau-lepton candidates, for the 
data recorded at ^/s = 161 — 172 GeV, the expected backgrounds, and a simulated Higgs 
boson signal with mh = m\ = 55 GeV. In those instances where more than one candidate 
pair passes the selection, the pair whose members have the lowest track multiplicity is 
chosen. If no distinction can be made, the candidate pair is chosen whose members have 
the highest value of the isolation parameter RlU'^cd- Here -Rem+cd is the ratio of the sum 
of the electromagnetic cluster energies and charged track momenta within a cone of 11° 
half-angle centered on the r candidate axis to that within a 30° half-angle cone. 

The hadronic part of the event, obtained by excluding the tracks and clusters from the 



selected r candidate pair, is then split into two jets using the Durham algorithm The 
invariant masses of the tau-lepton pair, m^,-, and of the hadron jet-pair, mhad? are calculated 
using only the tau-lepton and jet directions and requiring energy and momentum conservation. 
The resolutions of the mass distributions, later used in the calculation of exclusion limits (see 
Section are determined from signal events by fitting a Gaussian distribution in an interval 
which excludes non-Gaussian tails, resulting in a typical error of 6 GeV on the reconstructed 
masses. 

The numbers of observed and expected events after each stage of the selection are given 
in Table | for y/s = 161 GeV and 170 — 172 GeV. The agreement between the data and the 
expected background within the limited statistics demonstrates the adequate modelling of the 
selection criteria. The detection efficiency for a Higgs boson signal with mh = = 55 GeV is 
also given. Five events survive the selection while the background is estimated to be 1.29 events 
at a/s = 161 GeV and 2.54 events at i/i = 170 — 172 GeV. Figure ^ shows the positions of 
the surviving events in the (mT-T-,mhad) plane superimposed on the expected background. The 
detection efficiencies for various values of m,-,- and mhad, where m^,- is the mass of the object 
(h° or A°) decaying into the r-pair and mhad is the mass of the object decaying into the jet-pair, 
are given in Fig. |^. 

The efficiencies are affected by the following uncertainties: Monte Carlo statistics, typi- 
cally 2.2%; uncertainty in the tau-lepton identification efficiency, 3.6%; uncertainties due to 
modelling of selection variables excluding the tau-lepton identification, 6.3%; uncertainties in 
the modelling of fragmentation and hadronisation, 2.4%; and uncertainty on the integrated 



luminosity, 0.6% [^. Taking these uncertainties as independent and adding them in quadra- 
ture results in a total systematic uncertainty of 7.9% (relative errors). The uncertainties due 
to the modelling of selection variables, including those used in tau-lepton identification, were 
estimated by displacing the cut values by an amount corresponding to the difference between 
the means of the data and background Monte Carlo distributions. Using the same techniques, 
the uncertainty in the number of expected background events was estimated to be 32%, domi- 
nated by the systematic uncertainty associated with the requirement on | cos a, ■ cos aj \ , which 
is steeply falling for background at the position of the cut. Because of this predominance of 
one single selection variable in the uncertainty on the number of expected background events, 
the background is not subtracted when computing limits. 



4.3 The channel h^AO^bbbbbb 

When 2mA < the decay h*^^A*^A° is allowed and may be dominant. In these cases the 
process e"'"e~^h'^A°^A°A°A° can have a large branching ratio in the final state bbbbbb. Due 
to the presence of six b quarks in the expected signature, the events are characterised by a large 
number of jets and a large charged track multiplicity. To reduce backgrounds b-tagging plays 
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Figure 4: Selection variables relevant for the h°A°— ^qqr+r" analysis, (a) The energy of the 
most energetic electromagnetic cluster scaled by y/s after cut (3). (b) The pairwise isolation 
parameter (see text) after cut (4). The y/s = 161 and 170 — 172 GeV data are added together. 
The points represent the data. The shaded histograms show the qq background and the open 
histograms show the four-fermion background, normalised to the integrated luminosity of the 
data. Two-photon processes are not included. The dashed histograms are simulated signals 
for mil = = 55 GeV, where the displayed production cross-sections have been chosen for 
visibility. The background simulations are normalised to the integrated luminosity of the data. 
Arrows indicate domains accepted by the selection. 
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Figure 5: Efficiencies in percent for the signal process h°A°^qqr+r^ in the (Vn^T-,mhadj plane 
at a/s = 161 and 172 GeV, where m^-j. and m^ad are the invariant masses of the tau-lepton pair 
and hadron jet-pair, respectively. 
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Table 4: Effect of the selection criteria on data, background (normalised to the integrated 
luminosity of the data) and signal simulation (mi^ = = 55 GeV) at the two centre-of- 
mass energies for the signal channel h°A°— i>r+r'qq. The quoted errors on the background are 
statistical. 
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Figure 6: The position in the (mrr ji^had) plane of the li''A°— >-r^r qq channel candidates, 
superimposed on that of the arbitrarily normalised expected Standard Model background, for 
= 161 and 170 - 172 GeV. 
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a crucial role. The main background is from (Z°/7)*— i>bbg(7) with hard gluon emission. At 
161 and 170 — 172 GeV four-fermion processes also result in a small contribution. Backgrounds 
from two-photon processes are reduced to a negligible level by the event selection. 

The initial event selection follows that for h'^A'^^bbbb described in Section |4.1|. The data 



samples from ^/s = 130 — 172 GeV are all used. The following requirements are made: 



[1) The events must qualify as hadronic final states as described in p9 . 



(2) Jets are reconstructed using the Durham PBI algorithm with ?/cut = 0.0015. Events having 
five or more jets are retained. 

(3) As in the h'^A^^bbbb analysis, the radiative process e"'"e~— >(Z°/7)*-^qq7 is largely elimi- 
nated by a requirement on the effective centre-of-mass energy, in this case y/s' > 110 GeV. 
The distribution of the number of jets for events with four or more jets after application 
of this requirement is shown in Fig. |^(a) for data and simulated background and signal. 

(4) The number of charged tracks for the signal process is quite large, but the backgrounds 
from (Z°/7)*— i>qq(7) and four-fermion processes have long tails extending to high mul- 
tiplicities (see Fig. |^(b)). Candidate events are required to have more than 35 charged 
tracks. 

(5) Three or more jets are required to show evidence for b quark flavour, using the b-tagging 
algorithms discussed in Section || (BTAGl pBf and BTAG2 [0). The secondary vertices 
used in these methods are in addition required to have at least two tracks each with 
two /iVTX r-0 hits assigned 0. The decay length significances S, ordered in decreasing 
significance, must be successively S" > 8, 4, 3 for BTAGl vertices and S" > 8, 5, 3 
for BTAG2 vertices. Events must pass either the BTAGl requirements or the BTAG2 
requirements. The distributions of the significance of the most significant and the third 
most significant vertex for the BTAGl algorithm are shown in Fig. |^(c) and (d) for events 
passing requirement (3). 

Distributions of the variables relevant for the selection are shown in Fig. |^(a)-(d) for the 
130 — 136 GeV, 161 GeV and 170 — 172 GeV data combined. Within the limited statistics, the 
agreement of the data with the Monte Carlo simulations is reasonable, except that it was found 
that the PYTHIA simulation of e"'"e~— >(Z°/7)*-^qq7 underestimates the number of events 
with five or more jets. This discrepancy arises only in the Z° radiative return peak. The loose 
requirement on was made to decrease the number of Z° radiative return events in order 
to reduce the effect of this uncertainty on the final selection, while maintaining a high signal 
efficiency. This discrepancy, which would be eliminated with a harder requirement on v^, 
vanishes subsequently with the application of the remainder of the selection. 

The numbers of events passing each requirement, compared with estimates from the back- 
ground simulations normalised to the integrated luminosities, are shown in Table |^. Also 
shown are the detection efficiencies for simulated samples of e^e~— ^h'^A^— >A°A°A°— >bbbbbb 
with rrih = 60 GeV and m\ = 30 GeV. No events pass the selection requirements for any of the 
three data samples, consistent with the background expectations of 0, 0.23, 0.37 events for the 
130 — 136, 161, 170 — 172 GeV samples, respectively. 

The systematic errors on the detection efficiencies for the signal are dominated by the 
statistics of the 500-event Monte Carlo samples, typically 7%. In addition to the statistical 
errors, there are errors due to the modelling of the selection variables. A Monte Carlo sample of 
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Figure 7: Selection variables relevant for the h^A^^bbbbbb analysis, (a) The number of recon- 
structed jets using the Durham algorithm with ?/cut = 0.0015 after multihadronic event selection 
and application of cut (3), for events with four or more jets, (b) The charged track multiplic- 
ity after cut (3). (c) The highest secondary vertex significance for the BTAGl algorithm (see 
Section |^ for definition) after cut (3). (d) The third-highest secondary vertex significance for 
the BTAGl algorithm after cut (3). The distributions are added for the centre-of-mass ener- 
gies 130 — 136 GeV, 161 GeV and 170 — 172 GeV. Data are indicated by points. The shaded 
histograms show the qq('7j background, and the open histograms show the four-fermion back- 
ground, normalised to the integrated luminosity of the data. Two-photon processes are not 
included. The dashed lines represent a simulated signal with = 60 GeV, mp^ = 30 GeV for 
■s/s = 171 GeV, where the displayed production cross-sections have been chosen for visibility. 
Arrows indicate domains accepted by the selection. 
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Table 5: Effect of tlie selection criteria on data, background (normalised to the integrated 
luminosity of the data) and signal simulation (rrih — 60 GeV, niA — 30 GeV) at the three centre- 
of-mass energies for the signal channel h°A°— >bbbbbb. The quoted errors on the background 
are statistical. 
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10,000 events with rrih = 60 GeV, mp^ = 30 GeV for ^/s = 171 GeV was used to study systematic 
effects. Tlie most significant of tfiese effects on the efficiencies is the b-tagging requirement, 
where the modelhng of the b fragmentation and hfetime and the vertex finding algorithms can 
introduce systematic effects. These effects are similar to those for the h°A°-^bbbb analysis. 
The distributions of the significances of the first, second and third most significant vertices agree 
between data and Monte Carlo within the limited statistics with those for the (Z''/7)*— i^qq(7) 
simulation for the selection variables of this analysis. The systematic errors on the efficiency 
due to the b-tagging requirement were found to be 1.7% due to the impact parameter resolution 
and fragmentation uncertainties and 3.5% due to the b-hadron decay multiplicity uncertainty. 
The systematic errors due to the jet reconstruction, the requirement on V^, and the track 
multiplicity are 1.2%, 1.7% and 7.2%, respectively. The total systematic error on the detection 
efficiency, not including Monte Carlo statistics, is thus 8.5%. Additional systematic errors 
on the predicted total numbers of events arise from the error on the integrated luminosity 
(0.6%) and the theoretical uncertainty on the cross-section (1%). 

The systematic errors on the background estimates similarly include contributions from the 
modelling of the distributions used in the event selection and from Monte Carlo statistics. As 
mentioned previously, the PYTHIA generator underestimates the number of events with five or 
more jets by approximately 23% for the overall sample. Also the high-end tail of the charged 
track multiplicity distributions could be mismodelled. The background estimate is also subject 
to the modelling of higher-order QCD processes. In addition, the tagging of the third most 
significant vertex as a b is subject to the misidentification of jets which do not result from 
b quarks, as in the h°A°— i>bbbb analysis. To estimate all these sources of systematic errors 
related to the modelling of the Standard Model physics, the backgrounds calculated using 
the HERWIG Monte Carlo were compared with those from PYTHIA. The HERWIG Monte 
Carlo describes the jet rates of the data much better. The differences between PYTHIA and 
HERWIG in the predicted numbers of background events passing all selection requirements are 
—0.15, —0.03, 0.03 events for ^/l = 130 — 136, 161, 170 — 172 GeV and are consistent with zero 
within one standard deviation. These are less than the statistical errors on the predictions from 
either Monte Carlo generator, which for PYTHIA are 100%, 17%, 15% for the 130 - 136, 161, 
170 — 172 GeV Monte Carlo samples, respectively. The systematic errors due to the b-tagging, 
jet reconstruction, V^, and track multiplicity requirements were estimated to be 3.2%, 5.5%, 
4.4%, and 13.5%, respectively, using the same methods as were used to calculate the systematic 
errors on the detection efficiencies. Thus the systematic error on the background is dominated 
by the statistical error and amounts to approximately 100%, 23%, 21% for the 130 — 136, 
161, 170 — 172 GeV Monte Carlo samples, respectively. Since for this channel the predicted 
backgrounds are very small, no background subtraction is applied. 



5 Statistical combination of individual search channels 

The searches for Higgs bosons performed by OPAL have not led to any significant signals. The 
negative results in individual search channels, based on data at various centre-of-mass energies, 
are statistically combined to increase the sensitivity. A new statistical method [0], based on 
"fractional event counting" , is used for that purpose. The method is used to test the predictions 
of specific models (e.g. the MSSM or Two Higgs Field Doublet Model) for specific parameter 
sets by comparing them to the experimental results. 

The method assigns a weight to each candidate event for a given hypothetical Higgs mass 
m (test mass). The sum of the weights for all candidates is related to a probability to be 
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consistent with a signal-plus-background hypothesis and a background-only hypothesis in an 



analytical manner [T^. The probabilities for signal-plus-background and for background-only 
are related in turn, using Bayesian statistics, to a confidence level. 

The weight, Wjj(m), for each candidate j of each search channel i at the test mass m is 
determined by the product of a channel scaling factor, Ci(m), and another factor, fi{m,mij), 
which is determined by the expected mass distribution at m evaluated at the candidate mass, 
rriij. 

The channel scaling factors q are determined by the signal-to-background ratio: 

dim) = 1 + ^ \ \ / , 
^ ' \ ^^"^(m) ■ Si{m) j 

The expected number of signal events in channel z, Sj(m), is calculated using the model pre- 
diction for the cross-section and branching ratio, the integrated luminosity of the data set to 
which the search is applied,^ and the signal detection efficiency. The total expected signal rate 
is s(m) = I]jSj(m). The function Si{m,mij) is the signal probability density. Its maximum 
for any possible rriij is Sf^'^^{m). The function Bi{m) is the expected differential background 
rate per GeV for test mass m. For channels where the mass is reconstructed, Sf^^^{m) is in- 
versely proportional to the mass resolution. For channels without mass reconstruction, the 
ratio Bi{m) / S^^^{m) is replaced by the total background rate. This implies that the c, are 
larger for channels where the mass is reconstructed. 
The factor fi{m, rriij) is given by the ratio 



fi{m,mij) 



Si{m, rriij) 



The overall event weight for each candidate is thus given by: 

Wijirri) = K ■ Ci{m) ■ fi{m, rriij). 

The factor is a normalisation constant that fixes the largest value of Wij{m) to unity. 

The sum of candidate weights over all channels w{m) = J2i,j Wij{m) is converted to a 
confidence level CL{m) for this sum to be more likely due to signal-plus-background than 



due to background-only. For example for CL = 0.05, the signal hypothesis is rejected at 95% 
confidence level. 

In some channels as indicated previously, the background is statistically subtracted from 
the data to enhance the search sensitivity. In these cases, the subtracted background is first 
conservatively reduced by its systematic error. In calculating the Cj the best estimate of the 
background is always taken. 

The uncertainty on the signal detection efficiency is accounted for using the method de- 



scribed in 



There is an overlap between some search channels in selected signal and background events. 
In the case of overlap in the signal, there is the problem of assigning the correct weight to a 
candidate. This situation did not not occur for the search results of this paper, because the 
efficiency of a given search channel for the topology of any other channel is taken to be zero. 
Overlap in the background is only relevant when the background is subtracted from the data 
in areas where the same candidate is found simultaneously in different search channels. This 
situation does not occur for the candidate events found in the search channels used here. 



^ Identical search channels using different data sets are considered as individual channels. 
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Figure 8: Upper limits at 95% CL on (as defined by Eq. ^) using all SM search channels 
and assuming the SM Higgs branching ratios for the S° (dashed line), and discarding all search 
channels that use b-tagging but assuming a hadronic branching ratio of the ^ of 100% (solid 
line). 



6 Model-independent and Two Doublet Model interpre- 
tations 

Model-independent limits are given for the cross-section for the generic processes e^e~— S'^Z° 
and e"'"e~-^ S°P°, where S° and P° denote scalar and pseudo-scalar neutral bosons, respectively. 
The limits are conveniently expressed in terms of scale factors, and c^, which relate the cross- 
sections of these generic processes to those of well-known SM cross-sections (c.f. Eqs. (|l|), (@)): 

crsz = 4^, (3) 

^SP = A (4) 
Figure ^ shows the 95% CL upper bound for as a function of the S° mass, obtained from: 

atSZ 

^2 _ ^V95 



where Ng^ is the 95% CL upper limit for the number of possible signal events in the data, e is the 
signal detection efficiency and C is the integrated luminosity. The sum runs over the different 
centre-of-mass energies considered. The dashed line is computed using all search channels and 
assumes SM Higgs branching ratios for the S°. The solid line is computed assuming 100% 
hadronic branching ratio for the and uses only search channels that do not employ b-tagging 
and is therefore more generally valid. Below m^o ~ 5 GeV, the direct search loses sensitivity 
rapidly and the limit for is determined from Fzo only, as described below. 

Figure ^ shows contours of 95% CL upper limits for in the S° and P° mass plane, for the 
processes e"'"e~— ^ S°P°^bbbb and qqT~^T~ , respectively. In both cases a 100% branching ratio 
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Figure 9: Upper limits at 95% CL for (see Eq. (^) for: (a) the S*^f*— s>bbbb search channel 
assuming the bb branching ratio for both S° and P° to be 100%, and (b) the S°P°^qqr+r~ 
search channel assuming a 100% branching ratio for this final state. The invariant masses of 
the tau-lepton pair and hadron jet-pair are denoted m^-^ and mhad; respectively. The search 
efEciency is zero outside this area surrounded by the dotted line. 



into the specific final state is assumed. The contours are obtained from: 



with Ng^ being the 95% CL upper hmit for the number of signal events in the data. The 
results obtained for bbbb (Fig. ^(a)) are symmetric with respect to interchanging S° and P'^ 
while those obtained for r"'"r~qq are not. For this reason, the results for r"'"r~qq (Fig. §(b)) 
are presented with the mass of the particle decaying into r+r~ along the abscissa and that of 
the particle decaying into qq along the ordinate. The irregularities of the iso-c^ contours are 
due to the presence of candidate events that affect Ng^ . 

If the decay of the Z° into a final state containing S° or P*^ is possible, the width Fz 
will be larger than when only the known Standard Model decays are possible. The excess 
width that is still possible when subtracting the predicted Standard Model width from the 
measured Fz value can be used to place upper limits on the cross-sections of Z° decays into 
final states with S° or P° bosons. The additional width of the Z" resonance from final states 
not classified as lepton pairs can be extracted from the measurement of the branching ratio 
BR(Z'^— i^^^) = F(Z''— >£^)/F(Z''-^all). The limits obtained from this equation are not sensitive 
to radiative corrections that equally affect all final states, because these cancel in the ratio. The 
current measurement is BR.{Z^—>ii) = 0.03366±0. 00006 [|^. From this value and the predicted 



SM value obtained from the difference between the measured and predicted Z" width is 
(-1.2 ±4.4(exp) ±1.7(QCD) ±1.8 (EW)) MeV, resulting in an excess width less than 7.1 MeV 
at 95% CL. The uncertainties are due to experimental errors, QCD corrections and electroweak 
corrections, respectively. This excess width of the Z° can still be affected by corrections, such as 
vertex corrections and oblique corrections, that are different for different final states. However, 
these corrections are expected to be small for the Two Higgs Field Doublet Model (2HDM). In 
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Figure 10: Regions excluded at 95% CL in the Type II 2HDM. The black region is excluded 
using constraints from Fz only. The dark grey region uses the direct searches for the SM 
Higgs in addition, but discarding the search channels that use b-tagging, assuming a hadronic 
branching ratio of the h" of 92%. The light grey region is excluded for tan (3 > 1 in the 2HDM, 
assuming SM Higgs branching ratios for and A° . 



the case of a very light Higgs (mh < 2m^) that remains invisible or that decays into photons, 
electrons, or muons, some of the final states Z*S° with the Z* decaying to leptons could still 
be classified as a lepton pair and the above limit does not apply. However, in this case the 



decay-mode-independent search limit of is applied. 

In the 2HDM the bosons S° and P° are identified with h° and A°, and the couplings and 
are identified with sin^(/3 — a) and cos^(/? — a), respectively. The assignment of the possible 
excess width in Fz to the process Z°^h°Z* yields an upper bound for which depends only 
on the mass of h°; the assignment to Z°— ^-h^A^ yields an upper bound for which depends 
on the masses of both h° and A'^. Combining these limits, the black region in the mass plane 
shown in Fig. ^ is excluded at 95% CL regardless of h° and A° decay modes. In the 2HDM, 
disregarding the possibility of h°— >A°A°, the most important final states of the decays of the h° 
and A° bosons are bb, cc and r+r^. The branching ratios depend on tan/5, but the hadronic 
branching fraction always exceeds 92% For tan/9>l the bb channel dominates while for 
tan/5 < 1 the cc contribution may become the largest. 

In Fig. [ly the excluded area in the (mh,mA) plane is shown when the limits on and 
are combined. Below the dotted line, where the h°— i>A°A° decay is kinematically allowed 
and competes with the h"— >ff decay, the smaller of the detection efficiencies is used. The 
excluded area is therefore valid regardless of the h°— i>A°A° branching ratio. The dark grey 
area is excluded at 95% CL when BR(h°— i^qq)>92%, generally valid in the 2HDM. The limit 
in the 2HDM for equal h° and A° masses is at = = 41.0 GeV.Q The light grey area 
is excluded when SM Higgs branching ratios are assumed for h° and A°. This assumption is 
valid in the 2HDM for tan/5 > 1. In that case the limit for equal h° and A'' masses is at 

^ Throughout this paper numerical mass hmits are quoted to 0.5 GeV precision. 
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mh = mA = 56.0 GeV. The holes in the exclusion of the light grey area at the edge of the dark 
grey excluded region are caused by the single candidate event in the h''A''— *bbbb search. 



7 Interpretations of the search results in the MSSM 

In its most general form, the MSSM has more than one hundred parameters. In this paper 
we consider a constrained MSSM, with only five free parameters in addition to those of the 
SM. The model assumes unification of the scalar-fermion masses (mo) at the grand unification 
(GUT) scale, and unification of the gaugino masses (which are parametrised using M2, the 
SU(2) gaugino mass term at the electroweak scale) and scalar-fermion tri-linear couplings [A) 
at the electroweak scale. These simplifications have practically no impact on the MSSM Higgs 
phenomenology. In particular, a common scalar-fermion mass and tri-linear coupling is justified 
since only the scalar top (t) sector gives important contributions to Higgs boson masses and 
couplings. 

Other free parameters of the model are the supersymmetric Higgs mass parameter /z, tan f3, 
and the mass of the CP-odd neutral Higgs boson, ttia- As mentioned previously, the top quark 
mass has a strong impact on mh. Therefore, it is also varied within reasonable bounds. 

In scanning the MSSM parameter space, values of the above parameters are input to the 
HZHA event generator [^, ^ which is supplemented with parts of the SUSYGEN [^] program. 



The HZHA program provides the masses and couplings of all Higgs bosons as well as those of 
the supersymmetric partners. It also calculates the cross-sections for e"'"e~— >h°Z° and h'^A'^ 
at each centre-of-mass energy, corrected for initial-state radiation. SUSYGEN produces scalar- 
fermion masses on the electroweak scale, starting from the same input parameters. 
For the above parameters, the following ranges are considered: 

• mo: to 1000 GeV. The masses of physical scalar-fermions are obtained in SUSYGEN 
by running mg from the GUT scale down to the electroweak scale using the relevant 
renormalisation group equations. 

• M2: to 2000 GeV. The U(l) and SU(3) gaugino mass terms. Mi and M3, are calculated 
from M2 using the ratios of the corresponding coupling constants. Mi : M2 : M3 = ai : 
a2 : 03. 

• A: — 2.5-mo to 2.5-mo. This range is chosen to include all possible t mixings. 

• fx: -1000 to 1000 GeV. 

• mj^: 5 to 160 GeV. Beyond this range, values in 5 GeV slices around mA = 250, 400, 
1000 and 2000 GeV are also explored. Masses below 5 GeV are not considered since in 
that domain the decays of A*^ are uncertain. As discussed in the previous section and 
shown in Fig. |10|, for mh < 43.0 GeV, small values of mA are excluded on more general 
grounds. Note that the parameter mA which is used as input to HZHA and SUSYGEN 
is the "running mass", while the mass of the CP-odd Higgs boson A*^, which is used to 
express the results, is the physical mass. 

• tan (3: 0.7 to 50. This spans the theoretically favoured range, 1 < tan/3 < rrit/mi,, but 
also includes values less than unity, which are not ruled out by theory. 

• mt'. 165, 175 and 185 GeV. The range of values includes approximately two standard 



deviations of the measured top quark mass |12 
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Most parameters are scanned by dividing their ranges into bins of variable size and choosing 
the values in each bin at random ("stratified scanning"). Exceptions are the parameter nit, for 
which three discrete values are used, and the values of rrip^ greater than 160 GeV, for which 
5 GeV bands around the values are used. 

The following parameter scans are considered in order of increasing complexity: 

(A) The parameters mo and M2 are fixed at 1 TeV, fi is fixed at —100 GeV and mt = 175 GeV. 
The parameter A is fixed to (minimal t mixing) or a/6 TeV (maximal t mixing.) The 
parameters rriA and tan/5 are varied as described above. This simple scan serves as a 
benchmark corresponding to the scans proposed in |]TU . 



(B) We also consider two scans which correspond to cases of minimal and maximal mixing in 
the t sector, inducing small and large corrections to the h*^ boson mass, respectively. In 
these scans the parameters mg, M2, m^, tan/3, and mt are varied independently, /i is 
constrained to — O.lmq and the parameter A is set to 

• A = for minimal mixing, and 

• A = ^/6 niQ for maximal mixing. 

The mass mq of the "left-up" scalar quark at the electroweak scale is uniquely determined 



in terms of mg and the other input parameters [H9 . 



(C) In the most general scan considered here, all parameters described previously are varied 
independently in the ranges that are listed above. 

The number of parameter sets considered in scan (A) is about 50,000, that of the scan (B) 
is approximately 1,000,000, and that of scan (C) is close to 6,000,000. 

Before comparing the above parameter sets to the data, these are subject to a selection 
on the basis of theoretical and experimental considerations. Only those sets are selected that 
provide a ti mass larger than the lightest neutralino mass. Additional experimental constraints 
are applied, which do not relate directly to the searches described here. A parameter set is 
rejected if the sum of the partial decay widths of Z°— >-h°Z* and Z*^— »-h°A° is more than 7.1 MeV 
(see the discussion of the constraint from Tz in Section Parameter sets giving rise to chargino 
or neutralino masses or stop masses |51[] excluded by OPAL searches at LEP2 are also 



discarded. 

It has been shown that large values of A and yU may give rise to non-zero vacuum expec- 



tation values for the t fields which break charge and colour symmetry. Criteria for charge- and 
colour-breaking (CCB) minima of the MSSM Lagrangian have been determined but these 



may be substantially modified if also "tunneling" from the CCB minimum to the electroweak 
minimum is taken into account. A calculation that includes the effect of tunneling is available 
for one specific scenario out of a number of distinct possibilities [0 , but a complete treatment 



of CCB criteria is not yet available to our knowledge. A simple approximate criterion to avoid 



CCB minima is 52 



A^ + 3fi^ <x {mi^^ + mi^^), (5) 

where mi_^ and mi_^ denote the left- and right-handed scalar top masses and x~3. For the 
specific calculation that includes tunneling this bound was shown to be modified to x~7.5 . 
For the general scan (C) of the MSSM parameter space results will be shown without applying 
CCB criteria, and for CCB criteria applied with x = 3 and x = 7.5. 
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The experimental exclusion limits, at the 95% CL, are presented below, separately for the 
scans (A), (B), and (C). In each scan the total predicted number of events from all search 
channels is calculated using cross-sections, branching ratios, luminosities and search efficiencies 
for all different MSSM parameter sets. From this expected signal prediction and the number of 
observed events a confidence level is calculated according to the prescription in Section ^. For 
scans (A) and (B) the 99% CL exclusion is also shown, to indicate the sensitivity to the chosen 
exclusion confidence level. The results are presented for each scan in four figures: (a) in the 
(mh, rTT-A) plane for tan > 1, (b) in the same plane for tan /3 > 0.7, (c) in the (mh, tan /?) plane, 
and (d) in the (mA, tan/3) plane. For scans (A) and (B) the lower limits for the minimal and 
maximal mixing cases differ by only small amounts and only the weaker of the two exclusion 
limits is given. The theoretically accessible area corresponds to the larger one, for maximal 
stop mixing. The theoretically inaccessible areas are shown in the figures as speckled. 

The results for scan (A) are shown in Fig. |lT]. The structure near = 60 GeV is caused 
by the candidate in the analysis for the LEPl leptonic channel at mh = 61.2 GeV. From 
Fig. |ll](a), for tan/5 > 1, 95% CL lower limits can be obtained for mh > 59.0 GeV and 
TTiA > 59.5 GeV. When the tan (3 range is enlarged to tan/3 > 0.7 (Fig. |TT|(b)), the lower limit 
of mh is not affected, but no lower limit on mA can be given. Figures 0(a) and (b) show the 
region mA > 5 GeV. For mA < 5 GeV and mh > 43 GeV the value of sin^(/5 — a) is always 
very close to unity in scan (A). Therefore, in this region only the Z^h*^ production process is 
relevant for the exclusion, and the limit depends neither on mA nor on the A° decay modes. 
Hence, the limit for mh at mA = 5 GeV is also valid for mA < 5 GeV. Figure |ll](c) shows 
the projection onto the (mh, tan/3) plane. Large values of tan/3 correspond to mh~mA. In 
the (mA, tan/3) projection of Fig. |ll|(d) the strongest limit for mA of 89.0 GeV at 95% CL is 
obtained at tan/3~1.4. However, at 99% CL this limit drops considerably. 

Figure |T2| shows the results for scan (B). Differences with respect to scan (A) are due to 
the possibility of having lower mi values. This leads in general to modified couplings and in 
particular, for some parameter sets, to a strongly enhanced branching ratio for h°— *gg. The 
wider range of mi in conjunction with mt = 185 GeV leads to larger theoretically accessible 
regions. Despite these modifications, many essential features such as the limit for mh at very 
high mA, the limits for mh and mA near mh = mA and the absolute lower limit on mh remain 
unchanged. The absolute lower limit on mA for scan (B) is mA > 23.0 GeV at 95% CL, but 
the range 26.5 < mA < 55.0 GeV is excluded at 95% CL. The small unexcluded "island" at 
60.5 < mh < 61.5 GeV and 23.0 < mA < 26.5 GeV is due to the simultaneous effect of a 
large branching ratio for h'^^A'^A^ and the LEPl leptonic channel candidate. Because the 
theoretically allowed region at low mA and high mh is increased in scan (B) with respect to 
that in scan (A), there is an additional unexcluded triangular region near mh = 68 GeV with 
a minimum value at mA = 55.0 GeV at 95% CL. 

For the results of scan (C), shown in Fig. 13, all exclusions are at 95% CL. For this scan 
the exclusion regions are shown for three cases: Without the application of any CCB criterion 
in black, for the CCB criterion of Eq. ^ with x = 7.5 in black and dark grey, and with x = 3 in 
black, dark grey and light grey. The situation changes dramatically with respect to scans (A) 
and (B) due to the appearance of unexcluded parameter sets with small values of mh. These 
points are characterised by a large mass difference mA — rn^ and small sin^(/3 — a), hence the 
h°Z° production is suppressed (cf. Eq. (|ID). The h°A° production is kinematically allowed at 
LEP2 energies in most of the affected region, but the cross-section is small, and the current 
integrated luminosity is not sufficient to exclude these parameter sets. The existence of such 



parameter sets has already been pointed out in |55 
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mj^(GeV) m^(GeV) 

Figure 11: The MSSM exclusion for scan (A) described in the text of Section |^. Excluded 
regions are shown for (a) the (m-^, mx) plane for tan/3 > 1, (b) the (m^, mp^) plane for 
tan/? > 0.7, (c) the ('mh, tan (3) plane, and (d) the (rriA, tan (3) plane. The black areas are 
excluded at 99% CL. The black plus the dark grey areas are excluded at 95% CL. The speckled 
areas in (a ) and (b ) are theoretically inaccessible. 
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mj^(GeV) m^(GeV) 

Figure 12: The MSSM exclusion for the scan (B) described in the text of Section 0. Excluded 
regions are shown for (a) the (nih, mx) plane for tan/3 > 1, (h) the (m^^, itia) plane for 
tan/3 > 0.7, (c) the (rrih, tan/3j plane, and (d) the (rriA, tan/3j plane. The black areas are 
excluded at 99% CL. The black plus the dark grey areas are excluded at 95% CL. The speckled 
areas in (a ) and (b ) are theoretically inaccessible. 
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mj^(GeV) m^(GeV) 



Figure 13: The MSSM exclusion for scan (C) described in the text of Section ^. Excluded 
regions are shown for (a) the (rrih, m\) plane for tan (3 > 1, (b) the (m^, nix) plane for 
tan/3 > 0.7, (c) the (mi,, tan/?j plane, and (d) the (niA, tan/3j plane. All exclusion limits are 
at 95% CL. The black areas are excluded without applying any CCB criterion. When the CCB 
criterion of Eq. ^ is applied with x = 7.5 the dark grey areas are excluded in addition. For 
this CCB criterion with x = 3, the black, dark grey and light grey areas are all excluded. The 
speckled areas in (a), (b) and (c) are theoretically inaccessible. 
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If no CCB criterion is applied (the black area only is excluded), no absolute lower limit 
for mh can be given (Fig. |l3|(a) and (b)). For the CCB criterion applied with x = 7.5, the 
unexcluded area near = 5 GeV is due to the relatively weak limit on sin^(/? — a) for this 
mass range, as shown in Fig. |[ For x = 3, the lower limit on rrih is 50.0 GeV at 95% CL 
for rriA > 5 GeV in Fig. |l^(a), and it can be seen in Fig. |T3|(c) that this limit corresponds to 
tan (3 = 50. Such high values of tan jS always result in m\^'^mp^ when the soft SUSY-breaking 
masses are high (> 1 TeV.) However, as can be seen from the figures, in a general scan large 
values of tan/5 are not excluded for large mass differences \mp^ — mh|, due to large higher-order 
corrections involving a low-mass stop. For ttia < 5 GeV the model-independent limit applies (cf. 
Fig. 0). This leads, for the x = 3 CCB criterion, to > 43.0 GeV at 95% CL. Figure 0(c) 
shows that there is no exclusion at all at 95% CL in the (mh, tan/?) plane for the general scan 
if no CCB criterion is applied. This is due to the fact that for any (mh, tan/5) combination the 
other parameters can be chosen to simultaneously generate a large mA and a small sin^(/5 — a). 
However, many of these specific parameter sets are discarded by applying a CCB criterion. In 
contrast. Fig. |l^(d), which shows the (mA,tan/?) projection, is only marginally different from 
Fig. |TT](d) and|T2|(d) of scans (A) and (B). The exclusion limits on rriA are only slightly affected. 

The results in this section suggest that the MSSM parameter bounds, and in particular the 
limit on mh, derived from the benchmark scan (A) and minimal/maximal mixing scan (B) are 
not valid when a more general scan (C) of the parameter space is performed. 

8 Summary 

Searches for neutral Higgs bosons presented here have not revealed any significant excess be- 
yond the background expectation from SM processes. Evidence for both the e"'"e~— i>h°Z° and 
e"'"e~— >h°A° production processes has been searched for, allowing also for the decay h°— i>A°A°, 
when kinematically possible. Limits on these processes have been placed in a model-independent 
manner and within the framework of the 2HDM and the MSSM. These new limits substantially 
improve those previously published by OPAL. 

In the model-independent approach, limits have been placed on the product of cross-section 
and branching ratio for h°Z° production assuming SM branching ratios and fully hadronic final 
states, and for h'^A'' production with bbbb and r+r~qq final states. 

In the 2HDM interpretation, limits have been placed in the (mh,mA) plane both for the 
case of any value of tan/5 and for that of tan/5 > 1. Along the h°-A° mass diagonal 95% 
CL lower limits are set at mh = TnA > 4L0 GeV, independently of the value of tan/5 and at 
fnh = mA > 56.0 GeV for tan/5 > 1. 

In the MSSM, three different scans over the model parameters have been performed. For the 
simplest, benchmark scan, in which all parameters except ttia and tan/5 are fixed, a lower limit 
at 95% CL on mh > 59.0 GeV and mA > 59.5 GeV is derived for tan/? > 1. For tan/5 > 0.7, 
the limit on mh remains at 59.0 GeV, but no limit for mA is obtained. 

For the MSSM parameter scan with minimal and maximal t mixing, where only the /i and 
A parameters are constrained and all other parameters are left free, the 95% CL excluded area 
has slightly less extent than for the benchmark scan. 

In the general scan of the MSSM parameter space, parameter sets arise that cannot be 
excluded at 95% CL. All of these correspond to small, but non-negligible, cross-sections for the 
process e+e~— ^h°A°, and will either be observed or excluded as the integrated luminosity of the 
data increases. A fraction of those parameter sets can be excluded if requirements are made to 
avoid charge- and colour-breaking vacua. 
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